Multifractality in a broad class of disordered systems 



O 

o 

Oh- 
<D 
CO 

(N 



Olaf Stenull 

Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104 

(Dated: February 2, 2008) 

We study multifractality in a broad class of disordered systems which includes, e.g., the diluted 
x-y model. Using renormalized field theory we analyze the scaling behavior of cumulant averaged 
dynamical variables (in case of the x-y model the angles specifying the directions of the spins) at the 
percolation threshold. Each of the cumulants has its own independent critical exponent, i.e., there 
are infinitely many critical exponents involved in the problem. Working out the connection to the 
random resistor network, we determine these multifractal exponents to two-loop order. Depending 
on the specifics of the Hamiltonian of each individual model, the amplitudes of the higher cumulants 
can vanish and in this case, effectively, only some of the multifractal exponents are required. 

PACS numbers: 64.60.Ak 
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The physics of critical phenomena is usually character- 
ized by a few critical exponents. In certain, e.g. nonlin- 
ear, systems, however, the scaling behavior can be much 
richer. It can be even so complex, that its exhaustive 
characterization requires infinitely many critical expo- 
nents yJ . This is the famous phenomenon of multifrac- 
tality Ql ■ Recently studied examples are as diverse as 
heartbeat 0, quantum gravity Q and percolation M 
type problems like random resistor networks (RRNs) |fj 
and self-avoiding walks on percolation clusters • 

In this note we study the multifractal properties of a 
broad class of diluted physical systems, viz. those systems 
that can be described by lattice models with a Hamilto- 
nian that is a sum of bond energies U depending only 
on the differences = i?j — dj of continuous dynamical 
variables i?j and $j on the bonds between nearest 

neighboring sites i and j, 



H 



(id) 



(la) 



Here, 7^ is a random variable that mimics disorder. It 
is assumed to take on the values 1 and with respec- 
tive probabilities p and 1 — p. We focus on systems that 
are macroscopically isotropic and hence the bonds are 
assumed to be undirected, U{d) = U(— 1?). Moreover, U 
is assumed to have a well defined minimum about which 
it can be expanded in a power series in $ (the locus of 
this minimum is used to define 1? = 0). Otherwise, U is 
arbitrary. For example, it may be periodic or not. Given 
these assumptions and conventions, the bond energy has 
a Taylor expansion of the form 



(lb) 



with a\ being strictly positive. Diverse physical systems 
can be described by this type of Hamiltonian. The sim- 
plest of these systems is perhaps the RRN, where ??j 
corresponds to the voltage Vi at site i and is defined 
on the interval [—00, 00 ]. U(V) = icrV 2 , with a being 
the bond conductance, is the electric power dissipated 
on an occupied bond. A whole family of systems that 



can be described by the Hamiltonian H is the family 
of systems with x-y symmetry, i.e., systems that are in- 
variant under the orthogonal group O2 of rotations in a 
two-dimensional plane and the isomorphic group U(l) of 
transformations of the phase of a complex number. The 
most intuitive example here is perhaps a ferromagnet in 
which the spins are confined by crystal fields to lay in a 
certain plane. Other systems exhibiting x-y symmetry 
include superconductors, superfluid helium, the smectic- 
C and the hexatic-B phase of liquid crystals and so on. 
In the x-y model $i becomes the angle tpi that specifies 
the orientation of the spin at site i and is defined on 
the interval [— 7r, it]. The bond energy is a 27r-periodic 
function, U((p) = — K cos(</?), with K being the exchange 
integral. 

Since we are dealing with diluted systems, we are fac- 
ing basically a percolation problem. If p is small, there 
are only finite clusters. If p exceeds a certain threshold 
value p c , on the other hand, there exists an infinite clus- 
ter. At the threshold, p = p c , the system undergoes an 
isotropic percolation (IP) transition. Hence, the order 
parameter P^ (the probability that any site belongs to 
an infinite cluster) and the correlation length £ (the av- 
erage diameter of a finite cluster) scale as P^ ~ (p— p^ 13 
and £ ~ \p — Pc\~ u , respectively, where j3 and v are the 
well known critical exponents of the IP universality class. 
Here, we are interested primarily in physical processes 
taking place on the clusters like electric conduction or 
the interaction of spins. We will see that the cumulants 



Cf{x,x') = {(#l x ,)X 



(2) 



are adequate and convenient observables to investigate 
the multifractality of such processes. (• • ■ ) c stands for the 
cumulants of the average (■ • ■ ) with respect to the Hamil- 
tonian Q, e. g., {i9 2 ) c = (t? 2 ), {i9 4 ) c = {i9 4 ) - 3(i9 2 ) 2 and 
so on. {■••}<£ denotes averaging over all configurations of 
the diluted lattice and the star indicates the constraint 
that x and x' must be connected. 

To our knowledge, the scaling behavior of the cumu- 
lants © is not known to date with 2 exceptions: (i) 
the RRN where one has conventional gap scaling be- 



cause H is harmonic, C$\x,x') 



c'|^ /u , with 



2 



<p being the resistance exponent known to second or- 
der in e expansion and (ii) the diluted x-y model where 

C^\x,x') ~ \x-x'\^ u and cf ] (x,x') ~ \x - x'\^l v 
with <j) c being a critical exponent associated with correc- 
tions to scaling that is known to first order in e expan- 
sion [Tol ITU . The purpose of this note is to reveal the 
scaling behavior of the G^ for all systems covered by the 
Hamiltonian for all I. Using renormalized field the- 
ory we will explore the intricate connection of the present 
problem to the renowned noisy RRN |8| to arbitrary or- 
der in perturbation theory. We will show that the ci 
scale at the percolation threshold as 



the replica space Fourier transform ip^(x.) — fgexp(—i\ ■ 
6*)$(x, 9) of the order parameter, where A is the replica 
variable conjugate to 9. A takes on values on a discrete 
D-dimensional torus. The quantity "0x( x ) ^ s designed so 



that its correlation function 



G x,x';A =(^(x^_r(x') 



H 



(7) 



where (•••}« indicates averaging with respect to the 
Hamiltonian QJ, provides convenient access to the cu- 

mulants Cjjp . Applying a standard cumulant expansion 
one finds 



C i J ) (x,x') = A l 



r >\il>l/v 



(3) 



G 



A = 



exp 



with the exponents tpi being identical to the noise expo- 
nents of the RRN [2j and with the Ai being amplitudes 
which depend on the specifics of U , i.e., Ai ~ a;. 

Field theoretic model-In order to apply field theory 
and renormalization group (RG) methods we need to con- 
dense the Hamiltonian into a field theoretic Hamilto- 

.(0 



U=i 



(-1)' 

(20! 



(8) 



where K t (X) = ^JAW] 21 is homo geneous polynomial 
in A of degree 21. Equation ijSJ shows that C$ can be 
calculated via taking the derivative with respect to K\ (A), 
or in other words, that G(x, x'; A) is a generating function 

man that is suitable for studying the C { J> . This can be for the Cjjp . This property will play an important role as 

we go along; it will alow us to extract the scaling behavior 
of the C$ from that of G(x, x'; A) which in turn can be 
calculated by using field theory and RG methods. 

Diagrammatic perturbation theory.-As usual, the cen- 
tral element of our RG analysis is a diagrammatic per- 
turbation calculation. Rs constituting elements are 
the three-leg vertex — g and the Gaussian propagator 
G(k,A){l - Sxfi}, where G(k, A) = (r + k 2 + wX 2 )- 1 
and where k is a momentum or wave vector conjugate to 
x. Due to the factor {1 — 6? g}, which enforces the con- 
straint V'o( x ) = stemming from Jg-$(x, 9) — 0, the prin- 
cipal propagator decomposes in a replica carrying part 
G(k, A) and a part G(k, X)d^ g not carrying replica vari- 
ables. Each principal diagram decomposes into a sum of 
replica carrying diagrams consisting of these two types of 
propagators. 

Note that none of the irrelevant vi appears in the prop- 
agators. This is important because treating the vi in the 
same way as the relevant couplings r and w would ruin 
our perturbation expansion, i.e., increasing orders in an 
expansion of the Feynman diagrams in terms of the vi 
lead to increasing superficial degrees of divergence. It is 
mandatory to truncate this expansion, or in other words, 
we should treat the vi by means of insertions of the com- 
posite field (or briefly operator) 



done by following the seminal work of Harris and Luben- 
sky (HL) on the RRN and the diluted x-y model ^} with 
the result 

n = j ' d d x jf{l$jjr(v,v?)* + |# 3 J, (4) 

where the Gaussian kernel is given by 

oo D 

K (V, V S ) = r A - wAg -J2 v iH ( V «<°>)' • ( 5 ) 

1=2 Q=l 

The order parameter field $(x, 9) lives on a continuous 
d-dimensional space with the coordinates x. It is subject 
to the constraint /g-$(x, 9) — 0. The variable 9 is a 
replicated analog of the dynamic variable d and lives on 
a Z?-dimensional torus |12|. The physical situation is 
recovered in the replica limit D — > 0. The parameter r is 
proportional to p c — p, i.e., it specifies the distance from 
the critical point, w is proportional to a\ and vi ~ a;. For 
vanishing vi the Hamiltonian reduces to the original 
field theoretic Hamiltonian of HL. The vi are dangerous 
irrelevant couplings as far as the are concerned. This 
can be seen by performing a scaling analysis in the replica 
variable 9 that leads to 

C { J\x, x'; r, to, {v k }) = w l fi(x, x';r, {v k /w k }) , (6) 

where /; is a scaling function. This shows that the vi 
exclusively appear in the irrelevant combination vi/w l . 

However, it turns out the leading contribution to gI 
vanishes upon setting vi to zero, i.e., information about 
the leading scaling behavior of is lost by omitting vi 
and this is why the vi are dangerous. 

Physical contents- -To fully appreciate the physical 
contents of the Hamiltonian J3J it is helpful to consider 



Oi(x) = j 



$(x, 



D 

E 

a=l 



V^,)'$(x, 



(9) 



For the following arguments it is useful to employ 
the so-called Schwinger parametrization, i.e., to rewrite 
the propagators by using the mathematical identity 
G(k, A) = Jq°° ds exp[— s(t + k 2 + w\ 2 )]. Let us consider 
a generic replica carrying Feynman diagram with succes- 
sive single insertions of Oi in each of its replica carrying 



propagators. The A-dependent part of such a diagram is 
of the form 

v i y si Y R i (X) ex p ( w Yl s ^f) ■ ( 10 ) 

i {ft} ^ j ' 

Here, the summation £{«} ^ s a summation over some 
complete set of independent loop replica variables {k}. 
The summations indexed by i and j are taken over all 
the replica carrying propagators. Xj = Xj ({/?}, A), where 
A denotes an external replica variable, is the total replica 
variable flowing through propagator j. The summation 
over {k} can be simplified by a completion of squares in 
the exponential and eventually approximated by an inte- 
gration. This integration is Gaussian and hence straight- 
forward. Taking the replica limit D — * and using that 
Ki(X) is a homogeneous polynomial of degree 21 one ob- 
tains 

^5>c 4 ({s}f ^(A) +■•■ , (11) 

i 

where Ci({s}) is a homogeneous function of the Schwinger 
parameters {s} of degree zero that depends exclusively 
on the topology of our generic diagram. 

The remaining steps of calculating our generic Feyn- 
man diagram consist of integrating out the loop momenta 
and the Schwinger parameters. These steps are entirely 
analogous to those well known from the field theory of 
IP and are skipped here for briefness. For background of 
the methods involved here, e.g, dimensional regulariza- 
tion and minimal subtraction involving Laurent expan- 
sions in e = 6 — d, we refer to Ref . 01 • 

Beyond these standard procedures there is one intri- 
cacy involved here that warrants further comment. The 
ellipsis in expression stands for various terms each 

of which contains a homogeneous polynomial in A. The 
polynomials of the omitted superficially divergent terms, 
however, all have a higher symmetry than if; (A). If an 
operator which depends on A via one of the more sym- 
metric polynomials is inserted into one of the IP Feynman 
diagrams it can generate all sorts of polynomials, or for 
that matter operators, but it can never generate Oi. This 
feature distinguishes Oi and makes Oi a master opera- 
tor 6] whereas the other operators are just slaves. All 
slaves must be taken into account in the renormalization 
process and one has, at least in principle, to deal with en- 
tire renormalization matrixes instead of simple renormal- 
ization factors. However, these renormalization matrices 
have a particular, simple structure. Due to this simple 



structure, the scaling exponent of a master operator such 
as Oi is completely determined by a single element of the 
renormalization matrix. Hence, for the practical purpose 
of calculating a masters scaling exponent, the slaves can 
be neglected. 

Comparison to the noisy random resistor network- 
The noisy RRN is a generalization of the RRN in which 
the conductances of occupied bonds are random vari- 
ables. Thus, there are two types of quenched disorder in 
this model, viz. the dilution and the randomness of the 
conductance of individual bonds. To treat the two types 
of disorder simultaneously, Park, Harris and Lubensky 
(PHL) H introduced a variant of the HL model in which 
the role of the D-fold replicated 9 is taken by a (D x E)- 

fold replicated voltage 9. The perturbation theory for 
the PHL model can be performed essentially by following 
the steps described above. The only noteworthy differ- 
ence is that the field theoretic operators £> RRN leading 
to the noise exponents describing the current distribu- 
tion on the network contain instead of Ki(X) the poly- 
nomials X RRN (A) = E|=l[£a=l( A(a,/3) ) 2 ]* . where ^ 
is the replica current conjugate to the replica voltage 9 . 
Though different, Ki(X) and Zf RRN (A) share two pivotal 
properties. First, they are of sufficiently low symmetry so 
that the corresponding operators Oi and C RRN are master 
operators. Second, both are homogeneous polynomials of 
degree 21. Thus, in both cases the perturbation theory 
leads to expression 11110 [of course with either K\ (A) or 

Kf" (A)] up to unimportant differences residing in the 
ellipsis, i.e., up to different slaves. Therefore, the two per- 
turbation theories lead to identical results as far as the 
scaling behavior of the master operators is concerned. 

Scaling behavior. -Having made this observation we can 
draw on the noisy RRN, in particular on Refs. @, for the 
remaining steps. Eventually we are led for the correlation 
functions at criticality to the scaling form 

G(x,x';X) =| X -x'|- 2/3/!/ {i?o+ 

oo 

+ Y B mKi (X) |x - x'|*/" + ■••}. (12) 
l=i 

In writing Eq. i|12fl we have used that K\(X) = A 2 and we 
have set w = vi. The B's are expansion coefficients. The 
multifractal exponents ipi are identical to the noise expo- 
nents of the RRN and hence they are known to second 
order in s 0, 



e 313 + Z {3327 + 8/ [1556 + / (2076 + 881/)]} - 672 (1 + Z) 2 (1 + 2 /) 2 H(2l) 2 
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FIG. 1: Dependence of the multifractal exponents ipi on / in 
three, four and five dimensions. 

where H(n) = Y^l=x V^- Note that ipi = (j> and ifa — 4>c- 
Figure[I]plots the dependence of the ipi on I for several di- 
mensions. Our main result (0) follows immediately from 
the scaling form 1)12(1 by taking the derivative with re- 
spect to Ki(X) evaluated at A = 0. 

Moments vs. cumulants. -Knowing the scaling behav- 
ior of the cumulants C$ one might wonder about the 
corresponding moments 

M ( i\ X ,x') = {{dl x ,)Y € . (14) 

If the bond energy U is harmonic, than one readily finds 
by virtue of the relation (z9 2i ) = {2l)\/ (2 l l\){d 2 ) l c that 

mW(x,x')~\x-x'\W», (15) 

i.e., the moments display conventional gap scaling. The 
situation is much more intricate if U is not harmonic 
because then the higher moments correspond to compli- 
cated sums of products of the cumulants. We cannot 



prove, but it is not implausible that 

(n<^): fc V-i---'F-" fcW % as) 
^ k =i > c 

Provided that this holds, one is led back to Eq. i|15|) 
for the leading behavior of the moments in the limit 
\x — x'\ — > oo because ipi is a strictly monotonically de- 
creasing function of I and hence X^fcLi n kipk < with 
1 = Efcli^fc- 

Concluding remarks.-ln summary, we have studied 
multifractality in broad class of systems which includes 
the RRN and the diluted x-y model. The number of 
critical exponents ipi required to describe the scaling be- 
havior of the cumulants defined in Eq. J2J corresponds 
to the number of terms required in a power series expan- 
sion of the bond energy U. In the RRN, U is harmonic 
and hence the cumulants (J2J) show no multifractality. In 
the diluted x-y model infinitely many terms are required 
and one has true multifractality in this case. Note that 
only the first few ipi differ significantly from their large I 
limit ipac — 1. Hence, systems where the bond energy is 
not harmonic but when expanded features several terms 
beyond harmonic order will be hard to distinguish exper- 
imentally from systems with true multifractality. One 
might say that these systems are effectively multifrac- 
tal. The scaling behavior of the moments corresponding 
to the cumulants J5J) remains a challenging open prob- 
lem. We hope that our work stimulates experiments or 
computer simulations to decide whether these moments 
inevitably display gap scaling or not. 
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